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FLIGHT SIMULATION OF THE 
LONGITUDINAL MOTIONS 
OF ANOTHER AIRCRAFT 


. SUMMARY 


Simulation of the control-fixed longitudinal dynamic character- 
istics of the A4D-2, a Douglas jet attack aircraft, was attempted 
with use of a variable-stability airplane, Princeton University's 
North American Navion, and with telemetry for data presentation. 
The variable-stability installation was built as a joint project 
for both longitudinal and lateral-directional studies, and although 
the attempt at simulation of the A4D-2 was not successful, it is 
hoped that suitable guide lines have been established for similar 
studies. As a secondary objective, the effects of feeding back 
the airplane responses u,&% , 9 and @ were investigated. 

The project commenced in October 1959 at Forrestal Research 
Center, Princeton University, Princeton, New Jersey, and was 


completed in May 1960. 





FLIGHT SIMULATION OF THE 
LONGITUDINAL MOTIONS 
OF ANOTHER AIRCRAFT 


; INTRODUCTION 

Variable-stability aircraft and their use in flight research have 
received increased emphasis with the advent of sonic and supersonic air- 
craft. To date, much work has been done along these lines, primarily by 
Cornell Aeronautical Laboratory, Inc., and the NACA (now NASA). For 
longitudinal studies in particular, several aircraft including a C-45, 
an F-94, a B-26, and a T-33, have been variously configured by Cornell 
Aeronautical Laboratory with electro-hydraulic systems used to modify 
elevator and other auxiliary control surface responses, in order that 
studies could be made of: 

1) aircraft natural modes that are essentially aperiodic and 

convergent; 

2) optimum handling characteristics; 

3) proposed handling characteristics of an aircraft in design stages; 

4) 4inertia effects; 

5) phugoid damping. 
Similarly, the NACA has devoted much effort towards such studies, also by 


use of variously configured variable-stability aircraft. Refs. 1 through 


8 apply. - 


ko 





In this report, simulation of the longitudinal response character- 
istics of the A4D-2, the Douglas attack jet aircraft, was attempted 
with a North American Navion configured as a variable-stability air- 
eraft. Only specific response characteristics of the phugoid and short 
period were studied, not pilot opinion or handling qualities, since 
stick-force gradients and the like were not provided. 

The variable-stability installation built in Princeton University's 
North American Navion was constructed around a 3-axes, stick-steering, 
attitude-and-rate-sensing auto-pilot. Electrical signals produced by 
various transducers were summed in the auto-pilot summing circuits, and 
the control surfaces were moved by electrical servos through the ship's 
control cable system in proportion to these signals. The variable- 
stability installation was used to alter the normal stability derivatives 
of the Navion and to introduce additional or psuedo-stability derivatives, 
thereby changing the aircraft dynamic response characteristics. 

The auto-pilot installation was a joint project in which both longi- 
tudinal and lateral-direetional characteristics were studied. The other 
half of the joint project, lateral-directional, is described in Ref. 9. 
All studies were performed at the Forrestal Research Center, Princeton 


University, during the academic year 1959-1960. 





THEORY AND ANALYSIS 

The normal procedure in the determination of aircraft dynamic motions 
is to solve the simultaneous, constant-coefficient, linear differential 
equations of motion that characterize the transient response of Tieearee 
craft to a disturbance from some equilibrium flight condition. This 
procedure ie@ outlined in Ref. 10. In the detemmination of the desiza 
requirements of the variable-stability installation, the longitudinal 
equations of motion of both the Navion and the A4D-2 were compared as 
outlined below. 
EQUATIONS OF MOTION 

The longitudinal (stick-fixed) equations of motion may be written 


as follows, in non-dimensional time with the angles in radians: 


C 
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Now let the forcing function, Se take the following non-dimensional 
form with angles in radians; 
be = K, § S + Ki a ae + Ku + K), © (4) 
where K) » Ka? KR,» and Ky, hereafter wile pew.-merced to ae gain 


constants. Equation (4) may also be written as: 
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Abe 

+ (GS2C Rs (5) 
Ae 

where §,, u, & , @ and 6 are understood to be perturbation 

quantities. 


Substituting on into the moment equation: 


(Cyt +(Cy + “Cmg DK + (Cys, - d=) © = - Cag, Ky by 


i 2 Koaq - Z 
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Zo write the moment equation in real time, p ea. Transposing 
@ea collecting terms: 
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Expressing all the equations in a determinent form in real time: 
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GAIN CONSTANTS 
solutions to these equations of motion were obtained by assuming 


the solution of the form: 


Upon substitution of these variables into the equations of motion, 
together with their first and second derivatives, three homogeneous 
algebraic equations resulted. Further, expansion of the resultant 
determinant gave the standard quartic equation in \ » impleitiyes 


terms of the gain constants, as shown: 
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In terms of modified transfer functions, the numerator terms were 











found to be: 
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Evaluation of the coefficients of the characteristic equation, in 
“4 


terms of the Navion stability derivatives as shown in Appendix A, yielded 


the following relations explicit in the gain constants: 


Cy, =1 (16) 
Cz = 5.30 + 15.75 Kp (17) 
Co = 9.88 + 15.74 K) + 32.27 Kp + 15.73 Ky (18) 
Cy = 408 + .468 K) + 1.284 Ky - 1.942 Ky + 32.25 Ky (19) 
5092355 + 1-051 . Seen K, + 1.285 Ky (20) 


Also as shown in Appendix A, like constants of the characteristic 
equations of the Navion and the A4D-2 were equated in the following 


fashion, and solutions for the gain constants were found: 
; | 
(Chdin-o= (Cl) Mop. nav. 
i 
6s = 
( iia (C3) nop .NAV. 


and so on, where the notation "Mod. Nav." refers to Modified Navion. The 


gain constants were found to be as follows: 
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SIGNIFICANCE OF GAIN CONSTANTS 
As may be seen in Equation (6), additional airplane moments were 
added to the longitudinal moment equation. For instance, the Cy G3u 


term is implicitly composed of the two terms Cn, + Cm Se Ke, where 


Cn os K is a moment increment added to the original moment oa Then 
the individual moments are as follows: 
“ty + mge 3 = Cm, + ASm, (7 
Cot Cm§e Kl = Cm +t ACny (22) 
Cag + Cm ge Kh = Cag + ACy, (23) 
c +C ee . C + AC (24) 
mae Me - "ae mae 


Fxamination of these equations shows the following: 


1) for the Navion, C, =0, so AC, amounts to a pseudo static 


oe 


stability derivative; 


th, 


2) the same applies to ACn, 3 


3) since aircraft responses to « and @ already exist, ACh 
and ACs amount to altered stability derivatives and changed 
response characteristics of the aircraft. 

Signs of the gain constants iiateeg that K) through K, Supply de- 
stabilizing moments, while Ky supplies a stabilizing moment. A study of 
the magnitudes of the gain constants shows the following: 


1) K,: for a positive 15° increase in angle of attack, the variable- 


O 
stability installation must apply .6 of up elevator; 





2) K,: for a positive pitch rate of 5°/sec., 1.2° up elevator; 

3) Kz: for an increase in speed of 12 mph, .3° down elevator; and 
4) Ky: for a positive 15° increase in pitch angle, .05° down elevator. 
As was shown by Equations (13) through (15), the forcing function 


Se) appears in the numerators, 


Cm & rm Ke one (which is effectively Ce Gre 


but the gain,constants do not so appear. This is clearly the result of 
making changes only to the moment equation, and none to the lift and drag 
equations. The above indicated that with provisions only for modifying 
moments, no simulation could be made of: 
1) amplitude of oscillations in u, «& , 8, and 6 during transient 
responses, and | 
2) steady-state values of the above for say an elevator step input 
forcing function. 
It was seen, however, that this would not limit simulation of the transient 


response characteristics damping and frequency. 





PROCEDURE AND EQUIPMENT 


ANALYTICAL STUDIES 

It was decided that a complete analytical study should precede the 
flight simulation, with particular emphasis on the magnitude of aéceleratiee 
normal to the aircraft's flight path, and some @mphasis on the effects of 
the Piessibi ity of simulating aircraft transfer functions of the A4D-2. 
In order to accomplish the above, a computer program utilizing the Goodyear 
Geda Analog Computer was organized. The dynamic response characteristics 
of the AXD-2 were studied first, and damping and Puoquese ofweach ef thie 
two longitudinal response modes, short period and phugoid, were determined 
from the computer and compared with those obtained by analytical methods. 

For comparison, response characteristics of the unmodified Navion were 
also determined, as outlined in Appendix B. Then by use of the computer 
the simulation of the A4D-2 was accomplished, and the effects of vari- 
ations of gain constants were studied. Normal accelerations, primarily 
during the short period, were calculated, also as outlined in Appendix B. 

In addition to the above, studies were made of the effects of the 
gain constants on aircraft response as each individual gain constant was 
fed back with increasing magnitude from a zero value. 
DESIGN OF THE AUTO-PILOT 

Simultaneously with the initial theoretical studies of the variable- 
stability installation, physical design considerations of the auto-pilot 
variable-stability combination were examined. The first step was a flight 
test proaeen to determine control cable forces encountered in the Navion 


during maneuvering flight. Flight test instrumentation used in these 


10. 





flight tests was kept simple to avoid an excessive delay in the construction 
of the auto-pilot. Test instruments included a hand-held spring force gage, 
a glass tube-and-spring accelerometer, and elevator, aileron, and rudder 
control position indicators, which consisted of markings on bulkheads of 

the aircraft close to the respective control and pointers attached to the 
cockpit contfols. Since auto-pilot construction had already begun on 
Navion N91566, Princeton University's second Navion, N5113K, which is 
essentially an identical airplane, was used for the flight tests. 

Flight tests consisted of one flight where measurements were taken of 
control stick force per g, stick force per velocity, and rudder and aileron 
control force per deflection of each. Control stick force per g was measured 
in steady turns, stick force per velocity was measured over a range of 
velocity below and above the trim speed, 122 m.p-h., in 10 m.p.h. increments, 
and rudder and aileron control force per deflection were measured in steady 
Sideslips. One test altitude was used, 3000 ft., and one configuration: 
landing gear and flaps up, canopy closed, cowl flaps closed, mixture rich, 
1900 rpm, and 22 in. of manifold pressure. All tests were performed as 
outlined in Refs. 13 and 14. 

With the Navion on the ground, the hand-held spring force gage was 
used to determine control cable forces encountered during the above flight 
test work. Although measurements of control cable forces were made diffi- 
cult by a cross-control coupling between the aileron and rudder cables, it 
was believed that these forces could be subtracted out with sufficient 
accuracy for design purposes. It was subsequently accomplished in this 


manner. 


fel 





Cable force determination was the linking step to the determination 
of required servo pulley gearing ratios. The two basic considerations in- 
volved in the gearing ratios were whether the given servo-pulley combination, 
with 160 in. lb. design torque, could provide the required cable forces, and 
whether the speed of response of the servos was such that the airplane-auto- 
pilot combination would be dynamically stable. The latter consideration is 
outlined in detail in Ref. 9, and the result of both considerations led to 
the design and construction of a larger servo pulley for the longitudinal 
axis for better speed of response. Existing lateral and directional servo 
pulleys were found to be satisfactory. 

With the help of detailed blueprints and some parts of a Lear, Inc. 
electric auto-pilot specifically designed for the Navion, location of 
separate variable-stability auto-pilot components was decided. With limited 
space available in the Navion, and with the existing arrangement of control 
cables, the logical location of the three-axis servo installation was the 
center of the bulkhead at the after end of the baggage compartment. De- 
sign and construction of the servo support mountings and brackets commenced 
early in November. Design and construction of the electrical and electronic 
components got underway about the same time, as outlined in Ref. 9. 

Angle of attack and sideslip angle sensors were procured from the Naval 
Air Test Center, Patuxent River. By use of a probe, previously designed for 
Similar measurements, the angle of attack and sideslip angle sensors were 
mounted in tandem, approximately one and one-half chord lengths ahead of 
mene teading edge of the left wing, at the tip. The dynamic and Svar 


pressure probe, also previously designed, was mounted on the right wing-tip, 


2 





also approximately one and one-half chord lengths ahead of the leading edge; 
the transducer used for pressure (airspeed) measurements was located just 
inside the right wing tip bow. 

An electrical stick and rudder was designed and installed as described 
in Ref. 9, along with the feed-back gain potentiometers and their mountings. 
The balance of electrical and electronic equipment including telemetry was 
designed and was installed in the baggage compartment; see Ref. 9. 

When the servo and pulley installation had been completed, elevator 
servo cables were fed aft from the servo installation and tied into the ship's 
elevator cables at a small angle, with use of cable tie-ins provided by Lear; 
aileron servo cables were fed forward through a pulley system to the then 
existing aileron follow-up cables, and rudder servo cables were fed aft 
through a single pulley to the ship's rudder cables. Tension on the ship's 
cables was set at 30 lbs., as per Navion maintenance specifications, and 
tension in servo cables was set at 25 lbs. 

Servo clutch slip levels were set such that the safety pilot, with use 
of the ship's controls, could, for example, over-ride the elevator servo of 
the auto-pilot with 16 lbs. of control force. In terms of torque limits 
for over-ride, slip levels were set as follows: elevator servo, 160 in. lbs.; 
aileron servo, 100 in. lbs.; rudder servo, 60 in. lbs. These levels were 
designed to permit control of the airplane in the cases of either inadvertent 
engaging of the auto-pilot, or failure of the auto-pilot disengaging system. 

Electrical and mechanical installations were mostly completed by the 
middle of February, and a complete static (ground) test of the servo system 


was inaugurated. The static check included: 


ates 





1) setting the servo limit switches; 

2) with full control deflection, checking the slack sides of 
the control cables for tension; 

3) measuring the control force required to overpower the 
auto-pilot; 

4. observing the structural effects on the servo and pulley 
installation, and surrounding mountings, during maximum 
PTorcertegec., . 

5) observing all cables for excessive cable stretch during 

rapid, large deflections of all controls; 

6) checking the complete aircraft control system for proper 
and safe operation; 

7) inspecting the aircraft structurally after completion of 
the above; and 

8) testing for total weight and aircraft center of gravity 
location. | 

Upon completion of the static tests and electrical trouble-shooting, 

at the end of March, the aircraft was test-flown to determine possible 
discrepancies of the aircraft and the auto-pilot installation. After the 
test flights, the Navion was given a somewhat less stoned but thorough 
static check. 

CALIBRATIONS 

Calibrations consisted essentially of determining first the number of 

degrees of control surface deflection per volt input to the servo, form 


given servo feed-back gain, and second the magnitude of input signal from 


Lh, 





each transducer for particular potentiometer settings, e.g., the number 
of volts input from the angle of attack vane at 5° angle of attack, with 
a potentiometer setting of 0.1. The four sensors of the longitudinal 
channel, angle of attack vane, pressure (airspeed ) transducer, pitch gyro 
and pitch-rate gyro, were calibrated in this fashion, as was the electrical 
stick. Aircraft response quantities which were to be telemetered, the four 
above plus elevator servo position (elevator position) were calibrated also; 
this amounted to measuring the telemeter d.c. output voltage for an a.c. 
transducer signal input. 
FLIGHT TEST PROGRAM 

The flight tést program consisted of three flights which were necessary 
to trouble-shoot communications, telemetry, and auto-pilot electronics, and 
one flight during which test data was taken. Time did not permit further 
testing. 


The configuration tested was as follows: landing gear and flaps up; 
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canopy closed, oil cooler doors open, 1850 r.p.m., and approximately 20 in. 
of manifold pressure. Test altitude was 6500 ft., and airspeed was 110 m.p-h. 
indicated (120 m.p.h. true). The forcing functions utilized in the tests 
consisted of approximately .05 rad. elevator steps or impulses; these were 
provided by the auto-pilot with deflection of the electric stick, which 
had been calibrated to produce approximately .05 rad. of elevator deflection 
for full aft deflection of the stick. Also, during phugoid tests, the wings 
were held level by use of the auto7pilot rudder controls. 

First, dynamic characteristics of the basic Navion were tested with 
the auto-pilot and the elevator step forcing function; with telemetry, the 


dynamic response characteristics were recorded on magnetic tape by an 
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Ampex 309-C tape recorder and on time history traces by a Sanborn 150 
recorder. Second, with an arbitrary gain potentiometer setting of .15 
for the u feed-back, and with the auto-pilot and a .05 rad. impulse 
forcing function, the short period and phugoid motions of the modified 
Navion were tested and recorded. In similar fashion, them , 8, and © 
feed-backs wére tested. After the completion of each test run, the auto- 
pilot was disengaged, the aircraft was flown back to test altitude and 
test airspeed, and the auto-pilot was re-engaged. 


For reasons outlined in later sections, the simulation of the A4D-2 


was not attempted. 


iG. 





RESULTS AND DISCUSSION 


GENERAL RESULTS 

The primary objective, that of in-flight simulation of the A4D-2 
with the Navion and its variable-stability installation, was not 
accomplished because of the effects outlined below. The secondary 
objective, which developed during the theoretical analysis of the simu- 
mapion, that of an investigation of the effects of variation of the 
feed-back quantities on airplane dynamic response, was accomplished 
through a theoretical analysis, but was not verified by flight test data. 
However, it is felt that the information contained herein concerning the 
variable-stability installation and the attempted simulation should be of 
help in establishing guide lines for similar studies. 
ANALYTICAL RESULTS 

As outlined in some detail in Appendix B, a comparison of the dynamic 
characteristics of the two aircraft, the A4D-2 and the Navion, was made. 
For ease of comparison, these characteristics are repeated here; factor- 


ization of the two characteristic equations gave, in seconds: 


Phugoid Short Pertca | 
AkD-2 34.9 39.5 
Navion 62.6 
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From the above, it may be seen that the variable-stability installation 

had to be designed to change the dynamic characteristics of the Navion 

in the following way: decrease the frequency of the phugoid and the short 
period modes, decrease the damping of the short period mode, and increase 
the damping of the phugoid mode. As outlined in the section on theoretical 
considerations’, such changes were provided for by the use of the feed-backs 
u,e ,@, and 9, and the related gain constants were found to be of particular 
magnitudes and algebraic signs. 

At this point it should be re-emphasized that the characteristic 
equation of the Navion was made exactly similar to that of the A¥D-2 by 
use of the feed-backs, but, in terms of transfer functions, the numerator 
Germs were in no way changed. The transfer functions of the modified 
Navion then consisted of a constant comprised of Navion stability deri- 
vatives, a fourth order denominator term peculiar to the AkD-2, and a first 
or second order numerator term peculiar to the Navion. 

As outlined in Appendix B, the effects of this kind of matching were 
investigated. In a comparison of steady-state values, i.e., the magnitudes 
of the perturbation quantities u,x ,6, and © as tw, the same amount of 
elevator deflection used as a step forcing function produced markedly 
different results. For example, a .05 rad. elevator step produced 
Us, = -246 in the A4D-2, and u,, = 6.15 in the modified Navion. As noted 
in the Appendix, the latter is well beyond the small perturbation theory 
employed in this investigation. Examination of the individual stability 
derivatives involved in the steady state transfer functions showed why the 
steady-state values differed so markedly between the two aircraft. Numer- 


ator terms causing the difference include C,, and C; , and a denominator 
, a 
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term, C;, further complicated the picture: For example, Cm: of the 
Navion is approximately four and one-half times the value of that of the 
AMD -2; ori of the Navion is about twice that of the A4D-2; and in the 
denominator, C, of the Navion is about half that of the A4D-2. When 
these parameters are multiplied together, as is done in the transfer 
function, the steady-state value of u of the Navion, considering an 
equal amount of elevator step input, should be approximately 20 times 
that of the A4D-2. It should be emphasized that the principal offender 
in the steady-state transfer functions is the derivative “eed In gener- 
al, the same comments apply to the other pereurbacied quantities, «,0 
and 9, to varying degree. The above results were verified by computer 
studies and are tabulated in Table [. 

As mentioned in EQUIPMENT AND PROCEDURE, computer studies were also 
made of the dynamic characteristics of the A4D-2, the Navion, and the 
modified Navion. Dynamic characteristics such as frequency and damping 
of the oscillatory modes peculiar to the "three" airplanes are also listed 
in Table I, and generally good agreement may be seen, except in the case 
of frequency and damping of the short periods; accuracy in the determination 
of short period characteristics from a time history trace is difficult at 
best. Fig. Bl shows the computer wiring diagram used in the computer simu- 
lation of the A4D-2, and it may be noted that it is a true representation 
of the type of feed-back system used in the Navion variable-stability 
installation. 

Results of these computer studies also included transient values of 


the dynamic responses u,x ,O and @, for both the A4D-2 and the modified 


Navion. These results are also tabulated in Table 1. It may be noted 


\ 
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that the initial amplitudes of the response quantities to an elevator step 
forcing function of magnitude .05 rad. appear reasonable for the A4D-2, in- 
cluding normal acceleration(g's), but such quantities for the modified 
Navion are not only markedly different, but are well out of the realm of 
the small perturbation theory, as was true in the case of steady-state 
values. For the above reason, emphasis was shifted to the use of an in- 
pulse forcing function; as may be seen in Table I, however, even a .05 rad. 
impulse forcing function produced large initial amplitudes in oscillations 
of the response quantities, including 3.2 to 4.0 g's during the short 
period mode. 

In order to determine the cause of such large differences between the 
two aircraft in amplitudes of oscillations of the response quantities ar. 
airplane transfer functions were again examined, with emphasis this time 
On transient response. This procedure is also outlined in Appendix==aaa 
amination of the transfer function 5s) showed that the gain of the trans- 
fer function is composed of, in the numerator, Cams and Cp» and in the 
denominator, a constant 2, h, alee Evaluation of the gain showed that 
Ome, and J were the main causes of the large difference in pete between 
the AkD-2 and the modified Navion, viz., the gain of the Navion was about 
four times that of the A4D-2. Examination of the transfer function 4, 8) 
showed that the gain of the modified Navion was approximately seven times 
that of the A4D-2, and to further complicate the picture, the damping ratio 
of the second order lead term of the numerator of the Navion transfer 
function is approximately one-third that of the A4D-2, and well below criti- 


cal damping. These factors appreciably altered amplitudes of response. 
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Mien it may be seen that “mg. “Eee altered Scere response character- 
istics of the modified Navion, and additionally, the terms h and $ took part 
in this alteration. Examination of these terms and others will now follow, 
and later sections will touch on the same subject. It may be noted that a 
reduction of the ragnitude of the elevator step or impulse will reduce the 
amplitude San 6. tions of the response quantities of the modified Navion, 
and such reduction is linear with decrease in size of the step or impulse. 
However, there is some lower limit to this type of reduction; [ras ieee 
vious section it may be seen that the .05 rad. elevator step, Tor exampice 
Should be reduced by a factor of about 20 to .14°. Examination of this from 
a practical standpoint, considering cable slop and resolution of the average 
auto-pilot, shows that this is the right direction but the wrong path. 

me indicates that if C, 55 could be re- 


duced, response characteristics would be improved; a solution in this ’digec rien 


Fxamination of the effects of C 


would entail a small, auxiliary surface to be used as a moment-changer vice 
the Navion elevator. This will be discussed in more detail in a later section. 
Réduction of the amplitudes of the oscillations could be effected by increas- 
ing the magnitudes of h and > ; this would indicate that the airplane to be 
Simulated should have approximately the same mass and the same moment of 
inertia eee about the Y-axis as does the simulator airplane. It 
may be shown, however, that with a © feed-back, moment of inertia character- 
istics may also be simulated. Then, in the general case, 7 should be made 
approximately equal to that of the simulated aircraft. 

In the last direction indicated above, that of using 4 light air- 


plane to simulate a light airplane, resolution of the auto-pilot again 


becomes part of the picture. This phase was also investigated, and as 
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was expected, in the case that the aircraft and simulator aircraft are 
much alike in terms of stability derivatives and mass characteristics, 
the gain constants, K,....K),, tend to become extremely small; then the 
required magnitudes of the gain constants would tax an auto-pilot which 
had resolution approaching infinity -- or characteristics of a pure 
analog computer. 

Because of the above, and other effects outlined in following 
sections, and the very limited time remaining for further investigation, 
emphasis was shifted to the secondary objective mentioned previously: 

a qualitative investigation of the effects of variation of the fecd=uacr 
quantities. This was accomplished by means of root locii plots and was 
verified y computer studies. Only feed-backs of the algebraic sign de- 
signed for the A4D-2 were studied. Details of the root locii plots may 
be seen in Appendix B, and the general effects of the feed-back quantities 
are summarized in Table II. General comments concerning these effects are 
as follows: the principal effect of feeding-back angle of attack is to 
cause the frequency of the short period mode to decrease dnd the mode 
tends to become a convergent, aperiodic motion; there is little effect on 
damping of the short period, and frequency and damping of the phugoid. 

Feeding back pitch-rate has the effect of altering dynamic character- 
istics of both the short period and phugoid modes; at very high gain, the 
frequency and damping of the short period mode decrease, and the mode 
tends to become aperiodic, convergent. Further, at low gain the damping 
and frequency of the phugoid increase, at some medium gain the damping de- 
creases and the phugoid becomes oscillatory and divergent, and for very 


high gain the phugoid becomes aperiodic and divergent. 
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Unlike angle of attack, the velocity feed-back has the principal 
effect on the phugoid. For small gain, the frequency of the phugoid de- 
creases, damping remains about constant, and the phugoid tends to become 
aperiodic, convergent. With another small increase in gain, the phugoid 
becomes aperiodic, divergent; further increase in gain results in greater 
ervergence. In contrast, the damping and frequency of the short period 
mode are in general only somewhat increased. 

The pitch angle feed-back, the only feed-back designed for stabili- 
Zing the motions of the modified Navion, also has a strong effect on the 
phugoid. At low gain, the frequency of the phugoid decreases and the 
damping increases until the motion becomes aperiodic, convergent. The 
short period is effected also, since the frequency increases and the 
damping decreases; this is again a small effect, however, as in the case 
of velocity Peed-back. 

Although limited time did not permit an investigation of the effects 
of the feed-backs taken simultaneously, ay the magnitudes designed for the 
A4D-2 the separate effects of the feed-backs may be qualitatively summa- 
rized as follows: velocity and pitch angle feed-tacks decreased the fre- 
quency of the phugoid of the Navion, thereby changing the frequency to 
that of the A4D-2 phugoid. The angle of attack feed-back was used to 
decrease the frequency and pitch-rate was used to decrease the damping 
of the short period. Finally, both pitch angle and pitch-rate were used 
to increase the damping of the phugoid. 

FLIGHT TEST RESULTS 
As noted in the section on EQUIPMENT AND PROCEDURE, the in-flight 


simulation of the A¥D-2 was not attempted; some of the reasons therefor 
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are outlined in the preceding section. Flight testing that was done con- 
sisted of individual variations of the feed-back gains, also as mentioned 
previously. Although the aircraft, the variable-stability installation, 
the telemetering station, and associated equipment, were thought to be in 
good working order at the time of the test flight, numerous difficulties 
led to completely unsatisfactory flight test data. 

Problems that arose with the Navion may be enumerated as follows: 
during flight testing of the phugoid mode, it was noticed that engine 
speed varied approximately 150 r.p.m. between high and low speeds of the 
phugoid. It was felt at that time that the added moment from change of 
thrust was causing an alteration of the dynamic characteristics of the 
phugoid. This was particularly noticeable in the testing of the effects 
of the feed-back u; for a small elevator impulse, with resultant small 
amplitudes of oscillations of attitude and airspeed, and with resultant 
small changes in engine r.p.m., the phugoid was a stable oscillatory 
motion. At the same feed-back gain setting, however, and with a larger 
elevator impulse, the phugoid became an unstable oscillation, during which 
time the engine r.p.m. changed considerably. 

Another problem encountered with the Navion was the difficulty in- 
volved in trimming the aircraft prior to each test run. It had been found 
on a previous test flight that the auto-pilot trimming circuit was entirely 
too sensitive and could not be used to accurately trim the aircraft longi- 
tudinally. Since too little time then remained for changes of trim 
sensitivity, 1t was decided that the aircraft would have to be re-trimmed 


. with the ship's elevator trimming system between each run of the test 
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flight with the auto-pilot disengaged. This was done during the test 
flight, but this method did not provide the cure: the Navion trimming 
system was also found to be very sensitive, and, in addition, manifold 
pressure required to fly the aircraft at the test speed in a level con- 
figuration varied by as much as four inches of mercury. The source of 
this additional enone was not found. This difficulty in assuming the 
trimmed configuration naturally affected the phugoid test results. 

As will be discussed later, the sensitivity of the pitch channel 
servo amplifier had to be reduced to prevent oscillation or cycling of 
the elevator at the flight test gain potentiometer settings (feed-back 
gains). This reduced the resolution of the elevator servo. The low reso- 
lution of the elevator servo was evident during the tests of the @ feed- 
back; at low amplitudes of oscillation, there appeared to be little 
response from the elevator servo. During larger amplitudes df oscil- 
lation, the elevator was being moved in what was felt by the test pilot 
to be rather large steps, as indicated by somewhat rough control of the 
aircraft by the auto-pilot and seemingly large step-type changes in 
attitude. 

Additionally, some turbulence was encountered at the test altitude, 
and attempts were made to avoid flying over freshly plowed fields, and 
the like. Also, communications proved to be unsatisfactory, despite the 
use of a communications relay station between the Navion and the telemetering 
station. Finally, considerable "drop-out" of the telemetered signal was 


encountered at the telemetering station. 
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The result of the above was unsatisfactory flight test data, and 
there was no time left for repeated testing. 
THE VARIABLE-STABILITY INSTALLATION 

The variable-stability installation was built around a three-channel 
auto-pilot, which in turn consists essentially of three servo amplifiers 
and three servo motors used to position the control surfaces. The electric 
stick is used to generate voltage signals, these signals are summed with 
voltage signals from aircraft attitude and rate sensors, and the resulting 
voltage causes the servo motors to position the control surfaces. The 
installation is pictured schematically in Fig. 1. Considering the longi- 
tudinal channel, the servo motor positions the elevator through a bridle 
fastened to the ship's elevator cables, as may be seen in the schematic 
drawing, Fig. 2. Also, since the various sensors and the electric stick 
are excited by a.c. transformers, the signals generated have both magni- 
tude and sense. 

The gain potentiometers provided are used to divide the sensor signal 
voltages, and the portions of the sensor signal voltages summed and fed to 
the amplifier are dictated by the gain pot settings. The gain pot Getuame 
in turn, are a function of required gain constants. The algebraic signs 
of the gain constants are determined by the polarity of the exciting trans- 
former, e.g., for the A4D-2 simulation, polarity such that down elevator 
results from an increase in u. 

in addition, a trim circuit was installed. The trim circuit pro- 
vides a voltage of magnitude equal to but of phase opposite to summed 


voltages from the sensors, thereby permitting auto-pilot engagement. 
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To illustrate the elevator-movement sequence, assume that the initial 
input is a voltage generated by stick deflection. The servo and bridle 
move the elevator until the servo follow-up circuit generates a signal 
equal in magnitude and opposite in phase to the stick input signal. As 
the aircraft responds to the elevator deflection, voltages are generated 
ifethne pitch gyro, pitch-rate gyro, airspeed transducer, and angle of 
attack vane. These voltages, after being summed in the pitch channel sun- 
ming circuit, as may be seen in Fig. 3, are fed through the servo ampiiiiee 
to the servo motor, which deflects the elevator. The elevator augments or 
damps the aircraft response, depending on the phase of the summed a.c. 
voltage fed to the servo motor. 

A detailed description of the variable-stability installation is con- 
jeined in Appendix C, including pictures of the installation, déscripczom. 
of the structural and electrical components, and graphs of the control cable 
hemees which were the basis of the structural design. 

wearety features of the auto-pilot consist of the following devices 
an auto-pilot disengage switch is provided on the left console of the cock- 
pit for the safety-pilot, and a disengage button is provided on the electric 
Seperm for the test pilot. In addition, an electrical cut-off deviceme 
mrevided. If a large error signal is received at the point of inpui. 
wmescervo amplifiers, an electrical device breaks the circuit to tm 
servos, thereby disengaging the auto-pilot. When this occurs, another de- 
vice disengages completely the gearing between the servo motor and the 
servo pulley; the pulley then rotates freely on bearings only. Finally, 


as noted previously, servo clutches may be caused to slip, i.e., the 
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safety-pilot may manually override the servo motors with the ship's 
controls. 

Structurally, the auto-pilot variable-stability installation was 
found to be satisfactory in all respects except for cable slop in the 
aileron system. Some bending of the aileron servo cable pulley brackets 
was discovered during the static tests; this bending resulted in aileron 
cable slack at very high control force conditions. These pulleys and 
pulley brackets were the original ones used in the Navion aileron follow- 
up system, and because of position, the brackets were not readily ac- 
cessible for structural stiffening. 

AUTO-PILOT LIMITATIONS 

As mentioned previously, the gain constants required for the simu- 
lation were found to be of such magnitude that, for example, the auto-pilot 
had to be designed to furnish about 6 of elevator 4eriectionan Se 
change in angle of attack. Considering excitation voltage on the servo 
motor (corresponding to about 4° of elevator travel per volt), for a ise 
change in angle of attack, the resulting voltage (from the angle of attack 
vane) which was fed into the summing circuit was about .12 volts. The 
other gain constants, in general, resulted in similar magnitudes of volt- 
ages, with the exception of the pitch angle gain constant, which was an 
order of magnitude smaller. In order to provide such magnitudes of 
voltages, with the available sensor excitation voltages ( + 7.5 to + 30 
volts), and in order to provide gain potentiometer settings that were 
reasonable (towards the middle of the pot indicator scales), 15,000 to 


500,000 ohm end-resistors had to be installed in the individual sensor 
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Sircuits in series with the gain povenvtiomecers. (For comparison, gain 
potentiometers and sensor potentiometers have resistances of approxi- 
mately 5000 ohms.) | 

In a chain type reaction, considerable quadrature was produced in 
the summing circuit, the trimming circuit could not "interpret" the phase 
shift of the voltage, with the result that the voltage in the summing 
Sareuit could not be zeroed, and the auto-pilot could not be engaged. 
Quadrature was forthwith adjusted, the auto-pilot was engaged, but with 
ey change of gain pot setting or sensor position, the electric cut-onr 
disengaged the auto-pilot. Quadrature, then, was found to be a function 
Srenot only pot settings, but sensor positions as well. in additd#onjguas 
indicated, the electric cut-out could not "distinguish" quadrature from 
a large error signal. 

To relieve the above situation, the emphasis had to be shifted to the 
Second objective, investigation of the effects of variations of single 
feed-backs. The end-resistors were removed, but some quadrature effects 
were still found to be present. It was found that, for example, as the 
angle of the angle of attack vane was changed, a small amount of phase 
shift still took place; the explanation was that with the change of re- 
sistance of the angle of attack pot (associated with movement of the vane) 
Came also a change in a.c. impedance, causing phase shift. 

Ho relieve this same situation, and to thereby permit flighnt-tesriae 
Sie cut-out circuit had to be adjusted to a much less sensitive lever? 
and the sensitivity of the servo amplifier had to be decreased to prevent 


cycling of the servo motor. The reduction of sensitivity of the servo 
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amplifier reduced the resolution of the servo motor. The result was a 
barely acceptable longitudinal auto-pilot, considering the small changes 
of aircraft response which occur, for instance, during the phugoid mode. 

To further complicate the picture, the only suitable pressure (air- 
speed) transducer that could be located during auto-pilot installation 
had a non-linear potentiometer. This could be an additional explanation 
for aircraft-auto-pilot behavior during flight testing. 

As mentioned previously, auto-pilot resolution was an additional 
problem. Initially, resolution of the auto-pilot amounted to approxi- 
mately + see i.e., the auto-pilot could position the elevator within 
one-half degree of the commanded deflection. This resolution was found 
to be a function of the excitation voltage on the servo motor, and the 
magnitude of the signal voltage required to close the servo relay circuit. 
To improve resolution, the excitation voltage on the servo motor was in- 
creased to the maximum available in the system, + 30 volts. In addition, 
the peck-size of the relay was decreased with use of the servo-amplifier 
throttling control, or in other words, the amount of time that the relay 
remained closed was decreased. This effectively increased the peck-rate 
of the relay. These two improvements increased the resolution to 
about + fh. In an attempt to further increase resolution, a high- 
frequency sinusoidal voltage was superimposed on the relay, or, the 
"dithering" process. The result was unsatisfactory in that the elevator 
cycled at some sub-harmonic of the dither frequency. Finally, the sensi- 
tivity of the servo amplifier had to be reduced to prevent cycling of 


the servo motor, which was believed to be induced by electrical noise 


30, 





in the system; this reduced the resolution of the auto-pilot, but an 
absolute value of the resolution was not thereafter obtained. However, 
considering the 05° of elevator required by the 9 feed-back in par- 
ticular, it may be seen that even the improved resolution is an order of 
magnitude less than that required. 

In summary, the resolution and the quadrature problems were additional 
indications that a small, auxiliary surface should be used as the moment- 
changer, instead of the Navion elevator. In such a case, voltages would 
have been much larger, and quadrature problems of such magnitude could 
have been avoided. Another approach would have been to design the auto- 
pilot variable-stability system for d.c. summing circuit operation rather 
Pea a.C. 

CALIBRATIONS 

Quadrature also affected calibrations. With the end-resistors in 
place in the sensor circuits, for accuracy, quadrature had to be constantly 
adjusted during calibration for each range of pot settings. In addition, 
the sensor output voltages, as measured in the summing circuit, were so 
small that “they” could not be measured accurately with the existing volt- 
meters; this necessitated the removal of the end-resistors, calibration 
with other much smaller end-resistors, and a proportionate scaling of 
the calibrated quantities. Final calibrations may be seen in Appendix C; 
they are approximate for no end-resistors in place in the sensor circuits. 
WEIGHT AND BALANCE LIMITATIONS 

According to the Navion Operator's Handbook, maximum permissable weight 


is 2760 lbs., and maximum C.G. position is 30.5% of M.A.C., gear down. With 
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loaded configuration, including two pilots, full fuel load, and no passen- 
gers, the variable-stability Navion weighs 2750 lbs., -~ with gear down, 
the C.G. is 29.5% of M.A.C. 
ADDITIONAL DISCUSSION 

In summary, it has been shown in the preceding sections and in the 
Appendices, that with the Navion variable-stability installation, only the 
Navion moment equation was altered, not the lift or the drag equations. 
Further, provisions were made only for modifying two stability derivatives 
of the moment equation,C, and C,__, and for the addition of two pseudo- 
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derivatives to the moment equation, C, and C, . It has been seen that 


hy “< 

with such changes of the Navion moment equation, the characteristic. 
equation of the AXD-2 could be simulated, as was verified by use of the 
analog computer. Such simulation meant that the longitudinal dynamic 
characteristics damping and frequency of both modes of motion could be 
simulated. As was shown later, however, steady-state values and ampli- 
tudes of oscillations of the airplane responses could not be simulated. 
This last, coupled with auto-pilot quadrature and resolution problems, 
led to the abandonment of the simulation. 

Although the Navion variable-stability installation has a number 
of problems f the longitudinal channel, t#his*does not detract tremepae 
value of variable-stability installations in general. The literature 
contains the results of many successful installations, and the real 
value of such an installation is the fact that artificial control per- 
mits changes in the stability characteristics of the test aircraft while 


the aircraft is airborne. Without artificial control, it is impossible 
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to modify existing stability derivatives of a given airplane without 
major and costly structural changes. At the same time, without arti- 
ficial control, it is impossible to add new and additional stability 
derivatives, other than those normally contained in the three equations 
of motion, since there is no choice of what the airplane can "sense." 

With the variable-stability auto-pilot, the outputs of many types 
of sensing devices can be used, and there are practically endless possi- 
bilities in this direction, according to the literature, including Ref. 15. 
For example, in Ref. 4, the investigators used « and «x feed-backs to modify 
only the characteristics of the short period mode of the test aircraft; 
these response feed-backs were used in order not to appreciably change 
ime pnugoid. Also, u and u feed-backs have been used to only modify 
the phugoid response. As has been discussed in the present report, « and 
u feed-backs actually do just what is indicated above, and a quick ex- 
amination of the root locii plots will show that x and u also have the 
effects described above. 

Another approach to the variable-stability installation theory might 
be changes of the other two equations of motion, lift and drag, although 
this is not believed to be as practical as the moment-changer version. 

Also as indicated in the literature (see Ref. 4), the small, aux- 
iliary surface as a moment-changer is a practical appendage of the variable- 
stability installation. This has been verified by the results indicated 
in the present report. Instead of using a very powerful tool such as the 
Navion elevator to modify aircraft responses, it has been suggested, for 


conjecture, that the elevator trim tab could be used as the moment-changer. 
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Unfortunately, in the case of the Navion installation, this would mean 
additional weight for an additional servo motor, and other auto-pilot 
components, and it has been seen that maximum weight and C.G. conditions 


presently exist. 
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CONCLUSIONS 


ie view of the previous, it is concluded that: 


1) 


Simulation of the longitudinal flight motions of other 
aircraft cannot be successfully accomplished with the 
present Navion variable-stability installation. 

Considering simulation in general, unless all three 
equations of motion can be modified, steady-state values 
of aircraft responses and amplitudes of oscillations of 
responses cannot be exactly simulated; however, character- 
istic equations and thus dynamic response characteristics 
such as frequency and damping of the aircraft modes can 
be simulated. 

a) Considering gains of the airplane transfer functions, 
and auto-pilot resolution, a small, auxiliary surface 
should be used as the moment-changer in such type of 
variable-stability installation. Also, the test air- 
craft should have room for weight growth. 

b) Considering gains and numerator terms of the airplane 
transfer functions, in general, the simulator aircrage 
should have approximately the same time constant, +, 


as does the aircraft which is to be simulated. 
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3) 


4) 


In the present configuration the Navion is not suitable for 
in-flight testing of the effects of airplane longitudinal 
response feed-backs. 

Concerning feed-back quantities, «x has a primary effect on 
the short period mode, u and © have primary effects on the 
phugoid, and @ has effect on both the short period and the 


phugoid modes. 
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RECOMMENDATION 


Considering present Navion weight, present C.G. location, 
longitudinal auto-pilot limitations, and engine (r.p.m.) 
characteristics, it is recommended that future investigations 
with the Navion variable-stability installation be limited to 


lateral-directional flight testing. 
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Model 
Length 

Wing Span 
Wing Area 
Gross Wt. 
eG. Position 


Power Plant 


Propeller 


Landing Gear 


TABLE ITI 


NAVION SPECIFICATIONS 
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North American Navion 
teed SL. 

3308Cert- 

184.2 ft.° 
2750-1Ds. 
29.5% M.A.C. (with Auto-pilot) 
Continental, six cylinder, 
horizontal-opposed, 205 HP 
Hartzell, variable pitch, - - 


constant speed 


Tricycle, rétractaebie 








BPS a 
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FIGURE 2 
CONTROL SYSTEM SCHEMATIC 
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wae. control cables 
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APPENDIX A 
THEORY AND ANALYSIS 
MODIFIED NAVION DERIVATIVES AND BQUATIONS 


Given conditions: Clean configuration, power on (1850 rpm) 


h = 6500 ft., Q = .001957 slugs/ft.? 
Mes = 120 mph = 176 ft./sec. 
From Ref. 11 : No = .384 (Neutral Point) 
2 = ,095/deg. = 5.45/rad. 
Cae = 7 .025/deg. = -1.435/rad. 
as 2 
Cmag = 7 .15/rad. (Ty) sie = 3055 slug ft. 
C = - .O7/rad. 
Mao. / 
Assumption : That the derivatives above, which are ‘for the 


Ryan Navion, may be applied also to the North 
American Navion; both airplanes are standard 
“plueprint copies" made by two different air- 


craft companies. , 


Calculations 
Cm Cm = CL (Xoq - No) = 5.45(.295 - .384) 
= - .485/rad. 
w/S W/S = 2750/184.2 = 14.93 lbs./ft.° 
_ W/s | 14.93 = 
les “Ly = —q > TyaCools7y Maree = 193 
See 2 (a3) 5.45 
De Dy = Whe = “Fr (6.08) .85 = .333/rad. 





Calculations : D 





SL 493)" i 
De Pit). = Op  roeen ee) t TF (ONO ie 
ple a 2750 a 
a: a @ SV ~ 32.17(-001957)(176) (104.2) ~ 1. 
: _ _ 05 .5 = 
“_ a OSC ~ .001957(184.2)5.7 © 1aGe 
vor. 
"YG Yea 7 “qigee 
= 3055 - 85.5(1.39)°s 2890 slug ft.> 
Be 2 _tTIy _ 2890 _ 2 
Ky Ky ary “gees 22 ee 
“ | 
2 im 2 «33.85% 
h h “ae G2 7 a7 Tze = +05 
2 Cy #0 
“Dere ‘Derr “Dy 
C é: = iC 
Lor Lore =o 


Characteristic equation, Modified Navion: Dividing all coefficients 
by -h ae 
G, At + Gq A> OMe ee Cn eCESlO 
4 3 2 a O 


: Coeee)| 


ES STLaSy - -15 G - -O7 - .05(2.725 + .ohon)| 
5-30 + 15.75 Ky 


where Cy = 


2g) | 
o3 hf mag Ga ma 





C 
- hf Lex 
x 2 





i hc hc ' 
= c (c + Cola” - LQ) +Stlere(Cp Cr) 
hr? les. dex mae aw 2 a 


c : 
ar ta Cag 2) oe On x oo ‘en. 


Z OTL BEET [-ouon { - .Of - .15G5 - .05(2.785) 


i — (.333 - 493) + 2.725( - .156)} - -4856, 


= 1.1351 


= 9.88 + 15.74K, + 32.27K, + 15.73K, 


(CLg: Cn 


1 | Cr Cc 
= eee : apc? * Cag G)) + ~p tert (Cua Go 


C C 
eT) = ea £ 
on eae “Mg g Ga) 2 G3op,, + KuCn § age" 2) 


da Go 
2 
= -] oll 2.725( -.15)Go - ese) (493) 4 -.15G, 
~05(2.46) e 


e C - 
=| Cyaan rer (2 19)¢3 = “S ( .333)G3- isp (.otonse. 705) 


li 


408 + .468K, + 1.284K, - 1.942K, + 3262's: 


i Cric C Gc, - SL cia +n fe ae 
art |S [ett 1 To 3 Wem s | Deff 5 


ie ) 
a eff 
Sate 6: - C 


_ -2493 : 
~ 2(.05)3-32 [#930 485) @, 2.725 94 


1.435K), [ : 493 
ri OOM )e 725 - Meee 333 - Nos) 
(.05)3-32 - 


= .355 + 1.051 K, - 5-81 Ky + 1.285 Kj, 


A 3 





A4D-2 DERIVATIVES AND EQUATIONS 
From Ref. 12, the A4D-2 derivatives are as follows: 


Given conditions: Clean configuration, power on 








h = 6500 ft. V, = Mach .2 = 218 ft./sec. 
From Ref. 12: Cy = 3.62/rad. 
Cs A 
Cn g - .3265/rad. 
Cue 7 .145/rad. 
: OCm 
ee OSC peaks = — 
: 4 &) 
ev 
Cc 
C, = - 3-263/rad. -|° = 
q gc 
(96) 
Cn, = 0 
2 2 
Iy = 17,600 slug". ', S = 260 £2. 
qe = ]0.G9ft., We 107000 lbae 
CG = 259 MAC, b 227.5 it. 
Computation: 
_ WS. 10 q 
aes “Ly = ~q * I/B(.001957)(218)"200 ~ bo 
Cy * .059, e * .572, A = 2.91 
min 
c a2 
Do Cy = Cy + (Ci) 
O min TrAg 


308)2 
Seso4 ee SET = So 
1 sear 2( .828)3.62 
C : G = “~ = ° : = 1,1 ad. 
DEK Dx TA, r(2-91).572 7/ra 


. a | 10 ,000/g _ 
JS i osv = apie ars = 2.81 sec. 


Ay 





Computation: 


6 : =C. 2 aoe {ayes = - .00961/rad. 

Manx ° Cua x EDK" SV f- 2(218)2.81 / 
oe 7 7. _eeeeeso.s) 

Cn | Cag = Cn, -sUF = DSoime.sL OF .0288/rad. 


2 2 _l,  _ 17,600 _ 2 
v4 v — “"Tgs000/ Weteg wae 
. _ m_ _ 10,000/g ; 
fe = OSE = 7001057(BEO)IO.B = 26°6 
2 
-e@ yy 2 (Genome 
js “Ke 6 0 «5G.6 (LOSS ee 


y 
Characteristic equation, A4D-2: Dividing all coefficients by -hT ; 
i Tt t 
Cyt A #03 A> 40) A@ +c) A +c) = 0, where 
Cc, =1 
oo | - n¢c Pee fa 
3 hr? D =“) do ae 


-1 


Carey = .01715(.19 + 1.81) - .00961 - cose = 1.508 


t =i C CAC 

— +2 L 2 _ Ypvl 
C5 hf E be (Co Cy) Sot t} + Cp(Cmg,.+ Crag) 
CL 


>: 5 emag + a 


-l 


(01715(2.01)" | ox7as {anna at -.828) -9(2.2n} 


+ ,19(- .00961 -.0288) + 1.81(- .0288) 15] 


1.536 


ie 





; C C 
ae aoe es = Cugg + Cuma ) + ab CC -Crag ce ci 


-1 


= TOlTis)(2.B1)> .19(2.82x-.0288 -.145) + ath (628)(-~.00961) 
+.0288( 319) 


.0968 
|. 70,24 -( .828)°(-.145) ot 
Cy = Tipu gai 7 sees eae 


GAIN CONSTANTS 


Hevetinge the constants of the two characteristic eageatiedus. 


(Ci) aup-o = (Cyan amv. * 2 

C, = 1.508 = 5.30 + 15.75K, 

oh = 1.536 = 9.88 + 15.74K, + 32.27K, + 15.73K, 

C, = .0968 = 408 + .M6BK, + 1.2B4K - 1.942K, + 32.25K, 
C, = -O464 = .355 + 1.051K, - 5.81K, + 1.285K, 


Solutions of these equations yielded: 
K, = - .O407, dimensionless 


Ko 
K, = .0465, rad. 


oa 240, Sec, 


.0033, dimensionless 


Ky 
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APPENDIX B 
ANALYTICAL STUDIES 
COMPARISON OF DYNAMIC CHARACTERISTICS 
Factorization of the characteristic equations yielded the following: 
AMD-2, and MOD. NAV.: 
Phugoia NE oh O35ih 4 Moonee 


- .01755 + .18i 


= 
c 
Hi 


T /2 = 39-5 sec. P= 34.9 sec. 


Short Period Xe + 1.4734 + 1.452 =0 
ee iuaion = soot 
7} /2 = 94 sec. P= 6.59 sec. 
NAVION: 
Phugoid d= + .0212 \ + .0366 = 0 
Ty /p = 62.6 sec. P= 32.9 sec. 
Short Period + 5.28 A+ 9.72 = 0 


Ty /2 = .26 sec. P= 3.78 sec. 


COMPARISON OF STEADY-STATE VALUES 


The steady-state denominator term common to each aircraft transfer 





T@nction: 
C 1/2 (Cy - Cy) 3 
Dery ( Dy L L/2 
AN = Cr, 
SS. a 
Clore 2 0 
O Sik O 
ae cy? Cm 
o_ 


Bl 





Steady-state numerator terms: 





Cy, 
O 1/2(Cp - Cy) s 
C 
L 
(BJs.s. = 0 —= 0 
i Omg be Cm O 
G 
=C én 2a 
mee © fF te 
“Dere ° “1/2 
Mie) s.c. = Cees O 0 
| O - Cin 5 Se O 
2 
“L Ca ze 
ae 2) 98 e 
coer 1/2(Cp, - Cz) O 
oy 
(No)s.s. = Clore 5 = 0 





ALD-2: 
Evaluation of steady-state values for .05 rad. be 


(down elevator) step forcing function yielded: 


Us.,, = + -2h6 | Computer: + 3255 | 
Xs.s. = - -113 rad. = - 6.45° [Computer: - .114 rad.] 
Bs. = - .069 rad. = - 3.95° | Computer : - .O70 rad. | 


Be 





MODIFIED NAVION: 


For the same forcing function as above (including feed-backs): 


Be 5. = Gelb | computer: 5.8 | 
Xs is = - 1.1 rae. [_ Computer: - 1.0 rad. | 
985.5, = - 1.35 rad. [ computer: - 1.3 rad. | 


COMPUTER (DYNAMIC STUDY) 


AkD-2: Computer equations: 


. op ik Cy, 
e cee —— = ma? ¢ ot + — cama 
Lift : it putsplCp - Cra + 528 = 0 
, CL Cy, ; 
Drag : K+ Puts xX -@ =0 


pe eee (GR Yo + TBR + (Tie 4. =o 
Damping and Period: 
Phugoid T /2 = 40.0 sec., P = 35.0 sec. 
Short Period T /o = .8 sec., P 4% 8.0 sec. 


Transient values: For .05 rad. be (up elevator) 


step input, amplitudes of initial 


oscillations are as follows: 


uz - 45 

O 

Asp = 5-5 
fe) 

& ohu aT © 
89 = 19.5° 


B3 





Transient normal acceleration: 


Win - 1) tS 
: V 
or (n -l)g =a = po and @= 5 
then 6 =e) 
V 
’ V 
} [Ore ae Zz +1, From the compuver, 
6 mee rad./sec. initial amplitude of oscillation 
, during the short period; 
ug) 219( .08 
Dn =927 0+ = a8) 1 = 1352 


MOD. NAV.: Computer equations: 
Lift: u + .0297u - .0592« + .18250 = 0 
Drag: &% + .3655u + 2.02” -98=0 
Moment: © + .730u + 5.32% - .64la + 1.038& + .05180 + 2.226 
- 3.780 -"r> Koo. = 0 
Damping and Period: 
Phugoid Ti /2 = 39.0 sec., P = 35.5 sec. 
Short Period T = .5sec., P= 4.8 sec. 
Transient values: For .05 Ha. $ (up elevator) step input, 


amplitudes of initial oscillations are 


as follows: 


u = -~10 
A = +83 
n = 6.8 g's 


BY 





For .05 rad. & (up elevator) impulse 


forcing function, initial oscillations: 


6. Sie {Function of feed-backs | 
uo = - .75 
O 
Asp = 16 
, ° 
hu ~ 
Mo) $93.2 “ana Oaewac 


COMPARISON OF TRANSFER FUNCTIONS 

The denominator term common to all of the airplane transfer functions 
(not including auto-pilot dynamics) may be written as: 

A(s) =- np (S - A )(S - Ad)(S - A3)(S -A),), where 
Ay, are the roots of the characteristic equation 


(A4D-2 or MOD.NAV.). Numerator terms may be written as follows: 


C C C 

pa * &m, i. + Ba | ee 
: : atc 

a“ 





2 


m =-C St Ss + Ds + 1/2(2) | 
>, an me. © - 


2 oF = 2 
N, = - Cag beh S- + Cp + 5 8) S inxe (CpCy, - Cen + Cy) 


Now examine a typical transfer function, say 5.(8): 


| ie a 
Cm Cy S+=—S—_ 
Be (S) = me” 2% etch, 


BD 





Mmyvelueation of this transfer function gives; 
1.941 Ls + 2.90 | 


MOD. NAV.: = [0S AL). S =A] 


hoofs + .465 ] 
AND-2 : = “T(S -4))----(S -y)] 


“x 
Now examine $8): 
e 


2 C Cr, J 
Cng [S + 2055) 8 HGF) 


This may also be written in servo-mechanism notation as follows; 


2 2 
ey _* ls + 2808 +p ] where 
be [(s = Ay) Oh Stas (Ss ->i,)] 3 


K is the gain, § is the damping ratio, and W is the natural frequency 
of this minimum phase term (or second order lead term). 


Evaluation of the transfer function gives: 


MOD. NAV. : [ (6 = jee (s-d,) | 
pipe: = “Ts ae (s --,) J 


S) 
And examine (Ss): 


onal Cra, 1 2 
cms, Is + $ (Cp + “La\s 2 ape (CoCr, = CrCp a Cr } 
S) a (es 
fe) “ht 

e 
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MOD. NAV.: 


It 
r——7 
V— 
tA 
! 
a 
[+ 
=” 
—~ 
t) 
{ 
~” 
A toe 
ww 
L—! 


AkD-2 : 


| (SoA; eaten eee svocemiger > Me) 


ROOT LOCUS STUDY OF FEED-BACKS 
In a study of the feed-back quantities used to simulate the AD-2, 
the characteristic determinant for the Navion must first be examined. 


With angles in radians, and in real time, the determinant is: 


% cas Z 
(Cp + F,) 1/2(Cp) - CL) _ 
Cr (“Les Tp) “Sp 


With the addition of a u feed-back im the moment equation, the 
moment equation takes the following form: 
u + (Ci ee © o _ 
a ome. ( Tx Maa Sp) sg (Crs ies) tgee> = Cus be 


Then the characteristic equation takes the following form: 
hy 
~ "i (A-Az)(A- As) A - A3)(A - Ay) 





e 


Cr C 
eu, S a a. | -_ 


2Cy SF 
where A; Sesheeceetate Ay are the roots of the characteristic equation. 
eeecn roots naturally are modified by the root whith is part of the seeond 


term of the equation. The equation may be further simplified: 


BY 





a( \- Aj) greene etee ere (A -Ay) _ 


Ke = 0 
3 , 
(ha tebe 
Cp t 
where a = - ae De» and is a positive number. Then; 
© 
. Cree, ~ (.493)(5.45)  _ : 
Ae =" BCBS = 299 





2Cp S | NAVION 
Also, the roots of the characteristic equation of the Navion are; 


~ .0106 + .1944 


\ohu 


A 


S.p. 


- 2.64 + 1.671 


f 


where "phu" indicates phugoid, and "s.p." indicates 
short period. 

Berer to the root locus diggmem Fig. Be. For be positive, on 
fom elevator for increasing u, examine the ob condition: it may be 
seen that damping and frequency of the short period mode are somewhat 
increased (period and time to damp to one-half amplitude, P and Ty /22 
are decreased), but the characteristics of the phugoid mode are altered 
appreciably. The phugoid oscillation rapidly becomes unstable, then 
Deeomes aperiodic and divergent. 

In like manner, feeding back &, 8, and ©, singly and in order; 


results in the addition of the following terms to the original Navion 


Gmearacteristic equation: 


BS 





A = - .O147 + .2584 


e O 
a =- ht /Cry 3; O condition. 
c 


: a Cre C 
0: + Cm 6 S a [a2 + X(Cy op —)) + a en (Cp, = ou 


A= 0, - 2:01, + el 
—_—_— nt Cy, ; o condition 
C CL 
a oa 
@ : + rene ieee (Cp+ Ty | + — (Cp - c.) | 
, = - 2.01, + .O4; 180° condition 


her the results of the root locus studies, refer to Fuge uea. 
meeenen 5. For a summary of results of Appendix B, refer to Tables 


feomea IT. 


BY 
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APPENDIX C 


VARIABLE-STABILITY INSTALLATION 


The selection of the Minneapolis-Honeywell, USAF TYPE E-le, Auto- 
pilot imposed some limitations on the force and speed of response avail- 
able for control surface movement. Flight test results, presented in 
Figs. Cl and C2, showed that the force requifements of the vehicle for the 
variable-stability installation, Princeton University's North American 
Navion, (Fig. C3), were well within the 160 in. lb. shaft-torque output of 
the servo motors. The speed and frequency response characteristics were 
studied by another group working on the installation and are discussed in 
Ref. 9. 

Blueprints of an auto-pilot installation designed for the Navion were 
obtained from Lear, Inc. From this basic design two brackets were made 
of 1/8", 2h-ST, aluminum sheet and bolted to the partial bulkhead at the 
after end of the baggage compartment. Across these vertical brackets 
a 1/8", 24-ST, aluminum plate was bolted such that the plate lay in a plane 
approximately 90° to the ship's control cables. To support the servo motors 
a four-sided box, about eleven inches square, was fabricated of 1/8 ok -ST, 
aluminum sheet and 3/4" by 1/8", 17-ST, aluminum angle extrusions. To each 
of three sides of the box a servo was bolted, the bolts extending through 
aluminum angle stiffeners on the inner faces of the box sides. The remain- 
ing side of the box was bolted to the plate and bracket combination. Ad- 
ditional stiffening was supplied by 1/16", 24-ST, aluminum sheeting across 


the top and bottom of the box and two small 1/8", 2k-ST, sheet metal fillets 


Ci 





bolted to the vertical brackets and the partial bulkhead as lateral supports. 
The entire servo mount was over-—designed to preclude structural vibrations; 
the completed unit is shown in Fig. Ch. 

With the servos mounted vertically the capstans, or motor-driven 
pulleys, werg nearly parallel with the ship's rudder and elevator cables. 

. 

The bridle cable from the elevator servo, made of 3/32" aircraft control 
cable, was fastened to the capstan and fed directly aft about five feet to 
join the elevator cables at a small angle. The 1/8" rudder bridle cable 
was wrapped around the capstan and fed directly aft on the right side; a 
two in. fiber pulley was mounted on the under side of the sheet metal box 
to aline the bridle cable with the left rudder cable. The aitercna a 
of 1/8" cable, also was wrapped around its capstan, and the free ends of 
the bridle were fed forward through existing cable access holes in the 
bulkhead; a two in. fiber pulley was necessary for alinement on the right 
side. The 1/8" bridle cable was then mated to the 3/32" cross-control 
cable, the latter having been disconnected from the rudder cables. Standard 
aircraft cable thimbles and turnbuckles were used in all cable connections. 
The rudder and elevator bridles were attached to the ship's cables using 
sheet metal c-clamps. The bridle cable was united with the clamp by 
a 3/16" bolt passing through a thimble and the clamp; the clamp was tight- 
ened around the control cable with three more bolts. After completion of 
the bridle and cable matings, turnbuckles were adjusted for 25# tension 
in all bridle cables and 30# tension in the ship's cables. Servo clutch 
slip levels were adjusted to allow overpowering of the auto-pilot with 


reasonable control forces. The electrical servo stops were set to prevent 
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full deflection of any control surface in response to auto-pilot signals. 
A schematic diagram of the bridle-control cable arrangement 18 shown 
in Fig. 2. 

The conventional yoke and rudder pedals were removed from the right 
side of the cockpit. Two spring-loaded rudder pedals were mounted on the 
floorboards; a small cable attached to the pedals was used to drive 
a 360° potentiometer. The electric stick, a spring-centered flight con- 
troller with adjustable fluid damping, was bolted to a pedestal for appro- 
priate height. The flight controller had no provision for Stiek foree 
inputs; it had a small break-out force that was not objectionable. In 
addition to the above, Fig. 5 shows the auto-pilot trim box as it was 
mounted adjacent to the flight controller. 

Other modifications to the cockpit included the various auto-pilot 
end telemetry power switches mounted to the left of the instrument panel 
and the feed-back gain potentiometers grouped between the two front seats 
(see Fig. C6). The gain potentiometers are wire-wound potentiometers of 
approximately 5000 ohms each. With these feed-back controls eight air- 
craft responses ( Q, d, o. e, ®,% and u ) may be fed into the auto-pilot 
tO vary the response, and hence the dynamic stability, of the aircrare. 
Control surface travel can be varied by either the stick and rudder pedal 
gain potentiometers or by the servo follow-up gain controls. 

The major portion of the variable-stability installation was located 
in the baggage compartment just aft of the rear seat. The top shelf of a 
table was devoted to terminal strips that gave accessibility to the wiring 


of the system. These strips made it possible to monitor all the sensor 
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and excitation voltages; also, the sense of control deflection, for example 
down elevator for increasing u, could be changed at the terminal board 
merely by interchanging the transformer leads exciting the signal generator. 
The power to excite the primaries of the four multiple-wound transformers came 
from an inverter which converted the aircraft 28 volt d.c. supply to 115 

volt, single phase, 400 cycle a.c. 

The electronic portions of the auto-pilot were placed on a L6Tee shelf 
along with two rate gyros and the pitch-roll attitude gyro; three more rate 
gyros were mounted on the floor of the baggage compartment (Fig. C7). Al- 
though only three rates of response were measured a separate gyro was re- 
quired for each channel using the rate signal. The directional coupler 
feed-back loop and the control stick and co-ordinated turn capabilities 
of the auto-pilot were eliminated. An electronic safety device was ad- 
justed to disengage the auto-pilot, by releasing the servo clutches | 
altowing the capstans to rotate freely, if an error voltage of approxi- 
mately three volts suddenly appeared as an input to any servo amplifier. 

A disengage switch on the control stick was wired also, along with another 
disengage switch, for the safety pilot, adjacent to the main power switch. 

To record in-flight data an ASCOP pulse width frequency modulating 
telemeter station was used. The 15 channel filter conversion unit and 
the rotary sampling switch and transmitter unit were mounted in the 
baggage compartment (Fig. C4). The filter conversion unit allowed sampling 
of fourteen signals, reserving one filter for reference voltage. The 


rotary sampling switch operated at 20 cps and could sample 43 quantities 


plus two zero readings; hence a telemetry terminal strip was made to 
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parallel switch channels so that each of the 14 filter outputs could be 
sampled three times per cycle. Additional data on the internal components 
of the auto-pilot may be found in Ref. 9. 

To accommodate the remaining attitude sensors two booms were mounted, 
one at each wing tip rib as shown in Figs. C8 and C9. The booms, con- 
structed for previous testing, placed the sensing devices about four teer 
in front of the leading edge, ahead of the major pressure disturbances 
around the wing. Dynamic and static pressure lines extend from the 
pressure probe to a differential pressure transducer mounted inside the 
wing. Two Giannini yaw vanes were orientated 90° apart on the left probe 
to serve as sideslip and angle of attack vanes; vane mutual interference 
was considered beneath the resolution of the system. Since two position 
potentiometers were actuated by each vane only one angle sensor was re- 
quired to supply signals to both the yaw and roll channels. 

To complete the installation it was necessary to calibrate the entire 
system, that is to determine sensor signal voltages, for various gain 
potentiometer settings, throughout the operating ranges of the various 
attitude and rate sensors. Also it was necessary to determine control 
surface deflection for given input voltages. Then, to enable reading of 
the telemeter traces, ratios of a.c. sensor signal voltage to d.c. telemeter 
signal voltage had to be established for each of the quantities to be 
recorded. Results of the longitudinal channel calibrations are shown 


in Figs. ClO and Cll. 
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